Anova test
sq_df = read_csv('squirrel_tidy.csv') %>%
janitor::clean_names()
long_model_final = lm(long ~ shift + age + activity + reaction + sounds, data = sq_df)
anova(long_model_final)
## Analysis of Variance Table
##
## Response: long
## Df Sum Sq Mean Sq F value Pr(>F)
## shift 1 0.000250 0.00025024 4.2452 0.039446 *
## age 1 0.000446 0.00044617 7.5690 0.005973 **
## activity 1 0.000301 0.00030085 5.1038 0.023944 *
## reaction 1 0.000423 0.00042295 7.1751 0.007432 **
## sounds 1 0.001139 0.00113910 19.3243 1.141e-05 ***
## Residuals 3017 0.177841 0.00005895
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
lat_model_final = lm(lat ~ primary_fur_color + reaction + sounds, data = sq_df)
anova(lat_model_final)
## Analysis of Variance Table
##
## Response: lat
## Df Sum Sq Mean Sq F value Pr(>F)
## primary_fur_color 1 0.001256 0.0012564 12.034 0.0005297 ***
## reaction 1 0.001134 0.0011339 10.861 0.0009936 ***
## sounds 1 0.002101 0.0021015 20.128 7.512e-06 ***
## Residuals 3019 0.315195 0.0001044
## ---
## Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1